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(S| ' Abstract The coronal magnetic field is the primary driver of solar dynamic 

^ I events. Linear and circular polarization signals of certain infrared coronal emis- 

Q^ sion lines contain information about the magnetic field, and to access this in- 

"^ I formation, either a forward or an inversion method must be used. We study 

three coronal magnetic configurations that are applicable to polar-crown fila- 
^s^ ', ment cavities by doing forward calculations to produce synthetic polarization 

data. We analyze these forward data to determine the distinguishing charac- 
teristics of each model. We conclude that it is possible to distinguish between 
cylindrical fiux ropes, spheromak fiux ropes, and sheared arcades using coronal 
polarization measurements. If one of these models is found to be consistent with 
F^ ■ observational measurements, it will mean positive identification of the magnetic 

morphology that surrounds certain quiescent filaments, which will lead to a 
Q ■ greater understanding of how they form and why they erupt. 

^ I Keywords: Corona, Quiet; Magnetic fields. Corona; Polarization; Prominences, 

Ci ■ Models 






In order to understand coronal evolution and be able to predict dynamic events 
such as coronal mass ejections (CMEs) and solar fiares, we need to measure the 
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coronal magnetic field. However, making this measurement is not a trivial task. 
Much of this difficulty is due to the optically thin nature of the coronal plasma, 
and the relatively weak intensities of coronal emission lines. Measuring the field 
will bring resolution to many long-standing debates about the corona, including 
the magnetic-field morphology surrounding prominences. 

Prominences can be extremely stable on the solar surface; some polar crown 
filaments survive for many rotations (Gibson et a/., 2006), but they are also 
known to erupt suddenly. When prominences are seen on the limb, and are 
aligned with an observers' line-of-sight (LOS), they are often seen to be em- 
bedded in coronal cavities. These cavities are typically depleted in density by 
a factor of about two relative to the surrounding streamer (Fuller and Gibson, 
2009; Schmit and Gibson, 2011). Cavities are the coronal manifestation of the 
magnetic system that also includes the filament channel and the prominence 
(Hudson et al.^ 1999; Gibson et al.^ 2006; Heinzel et al.^ 2008; Gibson et al.^ 
2010). Because the cavity comprises the bulk of the system volume, the magnetic 
structure of the system can be determined from measurements of the cavity. The 
cavity-prominence structure is known to erupt bodily as a CME (Maricic et a/., 
2004; Regnier, Walsh, and Alexander, 2011). The initiation of these eruptions 
depends critically on the magnetic field threading through, and around, the 
prominence. 

Flux-rope models, with magnetic field wrapped around a distinct axis, and 
sheared- arcade models, without such an axis, have both been posited as possible 
morphologies of prominences and their surrounding magnetic field (see Mackay 
et al.^ 2010 and references therein). Flux-rope systems have also been explicitly 
compared to cavities (Low and Hundhausen, 1995; Gibson et al.^ 2006; Dove 
et a/., 2011). Both morphological models contain dipped field lines where mass 
can cool and condense into prominence material. Coronal polarization mea- 
surements could provide a means of distinguishing between these structures, 
which would lead to a greater understanding of not only the quiescent nature of 
prominences and cavities, but how they are formed, and how they destabilize and 
erupt. Our objective in this article is to determine the characteristic polarization 
signatures of these different models of prominence cavities. 

In Section 2 we discuss the difficulties of measuring the coronal magnetic 
field and the specifics of the Stokes vector in the Fe xill 1074.7 nm coronal 
emission line. Section 3 describes the forward calculations, followed by details 
of the three individual coronal models in Section 4. We analyze these synthetic 
forward-modeled observations and look for distinguishing features, which can 
be compared to solar observations of prominence-cavity systems to identify 
magnetic morphologies. We present these features in Section 5 and conclude 
with a discussion of our results in Section 6. 



2. Measuring the Coronal Magnetic Field 

There are several methods currently employed to determine the magnetic field in 
the corona. Given that the thermal conductivity along the field is high compared 
to across the field, the bright loops seen in extreme ultra-violet (EUV) and X- 
ray images of the corona trace out magnetic field lines. Using EUV images, it 
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is possible to follow these lines subject to projection effects (Aschwanden et al.^ 
1999). However, this method does not produce a measure of the magnitude of 
the field, and full three-dimensional (3D) traces require tomographic inversions 
and/or stereoscopic methods which can be tricky in the optically thin coronal 
plasma. Also, this technique only provides magnetic morphology information on 
specific bright loops, and not the full volume containing the magnetic field. 

Measuring Faraday rotation along a LOS to a known radio source supplies 
information about the LOS magnetic field, if the plasma density is known. This 
technique has been used to study the corona (Patzold et a/., 1987; Jensen, 2007), 
but a limited number of sight lines exist along which this technique is valid. 
Gyroresonant emission in radio wavelengths is related to the total magnetic 
field strength in the emitting region, but is limited to areas of strong magnetic 
field (> 200 Gauss) such as active regions (White and Kundu, 1997), although 
instrumentation with a broader spectral range, such as the Frequency Agile Solar 
Radiotelescope (FASR: Bastian, 2005), would allow for more extensive applica- 
tion. Observation of modified bremsstrahlung emission in radio and microwave 
frequencies provides information on the LOS magnetic field in the corona on 
the disk. These measurements sample a thin layer in the lower corona/upper 
chromosphere, not in the full coronal volume (Gelfreikh, 1994; Grebinskij et a/., 
2000). In some coronal emission lines, particularly in forbidden magnetic-dipole 
emission lines in the infrared, resonant scattering of anisotropic light in a mag- 
netized plasma can produce polarized emission subject to the Hanle and the 
Zeeman effects. We will now discuss the use of coronal emission-line polarization 
in more detail and further proceed to forward modeling of this emission. 

2.1. Coronal Stokes Vector 

Charvin (1965) was one of the first to show that linear-polarization signals 
from forbidden coronal emission-line transitions could be used to determine 
the plane-of-sky (POS) magnetic-field direction. Harvey (1969) was the first to 
attempt to use circular polarization to measure the LOS magnetic-field strength. 
Compared to modern technology, early detectors of the coronal Stokes vector 
had significantly lower signal-to-noise ratios, coarser resolution, and required 
longer integration times, which in turn limited temporal resolution. One such 
early detector was the Coronal Emission Line Polarimeter (Querfeld, 1977), 
which was a scanning photodiode polarimeter. A full- Sun measurement with 
this instrument would typically take about two hours and contain 1408 data 
points from 1.01 R© to 1.65 R© (Querfeld, 1977; Arnaud and Newkirk, 1987). 

In the intervening 40 years, there has been steady progress in the field of 
coronal polarization. Today, there are two main coronal polarimeters currently 
in use: The first is the Optical Fiber-bundle Imaging Spectropolarimeter (OFIS) 
on the Solar Observatory for Limb Active Regions and Coronae (SOLARC) at 
Mt. Haleakala (Lin, Kuhn, and Coulter, 2004). It generates 128 spectra from 
a 16 X 8 fiber optic array that subtends 5 x 2.5 arc minutes. It is capable of 
measuring full Stokes profiles at each of the 128 positions. This instrument has 
been used to successfully measure the linear-polarization strength and direction 
in the Fe xill 1074.7 nm line, and to determine a LOS field strength from 
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circular-polarization signals above an active region (Lin, Kuhn, and Coulter, 
2004; Liu and Lin, 2008). The second instrument is the Coronal Multi-channel 
Polarimeter (CoMP: Tomczyk et al.^ 2008), which is installed at the Mauna Loa 
Solar Observatory, and began taking full-corona measurements in October 2010. 
CoMP is an imaging coronagraph polarimeter with a tunable birefringent filter 
capable of detecting the Fe xill 1074.7 nm and 1079.8 nm lines as well as the He I 
1083 nm line. The new CoMP observations provide, for the first time, daily full- 
Sun observations of the magnetic field in the corona. The primary observables 
of CoMP are the four Stokes parameters (/, Q, /7, V). 

These observations are taken above the solar limb in the corona, which is 
optically thin in these wavelengths, and thus the measurements contain infor- 
mation from and extended LOS source. The polarization signal strength is weak 
compared to the line intensity (linear polarization/intensity ^ 10~^ and circular 
polarization/intensity ~ 10~^ for a one-Gauss field: Arnaud and Newkirk, 1987; 
Lin, Penn, and Tomczyk, 2000). It takes on the order of a few minutes to obtain 
a useable full- Sun linear-polarization measurement with CoMP, and circular- 
polarization measurements of sufficient signal-to-noise are made by averaging 
over an hour of data (Tomczyk et al.^ 2008). The polarization is the result of 
resonant scattering of anisotropic incident radiation by highly ionized coronal 
plasma in the presence of an external magnetic field. However, different aspects 
of this unified process dominate the linear- and circular-polarization signals of 
the coronal emission lines {e.g. Casini and Judge, 1999; Casini, 2002; Judge, 
2007; Rachmeler, Casini, and Gibson, 2012, and references therein). We restrict 
the rest of this discussion to the Fe xill 1074.7 nm coronal emission line. This is 
also the line used in the forward calculations. 

The linear-polarization signal is completely dominated by the Hanle effect: a 
depolarization of scattered light associated with a radiation-induced population 
imbalance of the atomic levels (Trujillo Bueno, 2001). The atomic alignment 
[a] describes this population imbalance. The transverse Zeeman effect, which is 
due to the energy splitting of the magnetic sublevels by the coronal field, is a 
secondary source of linear polarization, because the field strength is small (the 
transverse Zeeman effect is quadratic in the field strength) compared to the 
thermal width of these coronal lines. The Larmor frequency is larger than the 
inverse lifetime of the excited state, so the linear polarization signal occurs in 
the strong- field regime, also known as the saturated Hanle effect. In this regime, 
the linear-polarization strength and direction is dependent on the angle of the 
magnetic field, but no information about the magnitude is contained in the 
signal. 

The strength of the total linear polarization [L = \/Q'^ + If^] (same as P 
in Dove et a/., 2011; Rachmeler, Casini, and Gibson, 2012), is dependent on 
the angle [9] between the LOS and the local magnetic-field vector. Specifically, 
L oc sin^ Q such that L is strong when the magnetic field is in the POS, and 
weak when the magnetic field is along the LOS. The relative strengths of Q 
and U are used to determine the azimuth angle [^; U/Q = tan 2^], the POS 
angle of of the LOS integrated magnetic field. There is a 90° ambiguity known 
as the Van Vleck effect (van Vleck, 1925; House, 1977) such that the magnetic- 
field direction could be parallel or perpendicular to the measured ^. When the 
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local magnetic field is at the Van Vleck angle of roughly 54.7° with respect 
to solar radial, the light becomes unpolarized, and the strength of L goes to 
zero. When the magnetic field is less than 54.7° from radial, ^ is parallel to 
the direction of the POS component of B, but switches to perpendicular when 
that angle is surpassed (see, e.g. Figure 5(c)). The Van Vleck effect results in 
linear polarization directions in the corona that are mostly radial (Arnaud and 
Newkirk, 1987). If the location of the Van Vleck inversion can be identified, the 
90° ambiguity can be removed, although a 180° ambiguity remains. 

A measure of the magnetic-field strength is not possible with linear polariza- 
tion in this regime, but the circular polarization does contain information about 
its magnitude along the LOS. The Stokes V profile is proportional to BcosQ. 
The longitudinal Zeeman effect (which is linear in the field strength) is the main 
contributor to the circular-polarization signal. However, the atomic alignment 
can yield a significant correction to this signal, changing the amplitude of the 
anti-symmetric V profile, and therefore affecting the diagnostics of the magnetic 
field strength (Casini and Judge, 1999). 

Stokes /, Q, /7, and V are all dependent on the plasma parameters in the 
emitting region. They are weakly dependent on the temperature as long as the 
emission line is excited. All of the Stokes components are directly weighted by 
the density. At a given location along the LOS, the density dependence cancels 
when analyzing the relative polarizations [L/I and V/I] but this is not the case 
in a signal that is integrated along the LOS. For an integrated measurement, 
the signal will be dominated by those areas along the LOS that have the high- 
est density. Since collision tend to equalize the sublevel populations, a density 
dependence also enters in to the Stokes vector through <j. 

There are two general methods for interpreting the coronal-polarization mea- 
surements. The first is by inverting the signals into physical properties of the 
magnetic field. This is the approach taken with photospheric polarization data. 
However, because the plasma is optically thin in the corona, the signal is coming 
from an elongated source along a LOS. Inversions generally solve for a single 
point of emission, so not all of the calculations will converge to a solution. Inver- 
sions of these polarization signals require numerous initial assumptions about the 
emitting plasma, and are quite difficult due to multiple integrals that must be 
inverted. These calculations are known to be ill-posed (Judge, 2007). Information 
about atomic level-populations, and hence the plasma parameters, at each point 
along the LOS is required to solve the POS magnetic-field direction. In order to 
determine the field strength and direction everywhere, tomographic inversions 
are needed. The tomographic inversion process requires multiple viewpoints of 
the field, and if only one is available, as is currently the case, solar rotation 
must be used to generate these viewpoints (Kramar, Inhester, and Solanki, 
2006; Kramar and Inhester, 2007). This adds the additional assumption that 
the coronal field does not change appreciably over rotational timescales. An 
alternate approach to extracting information from coronal-polarization data is 
forward modehng. 

Our forward technique involves creating simulated polarimetric observables 
from models of the corona (Judge and Casini, 2001; Judge, Low, and Casini, 
2006). In the work presented here, we use this technique to study the differences 
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between several pre-CME magnetic morphologies, and expand upon the work 
begun by Judge, Low, and Casini (2006). The ultimate goal of this research is to 
determine if it is possible to use coronal polarization to positively identify flux 
ropes, or other magnetic morphologies, in the cavities that surround pre-CME 
filaments. 



3. Description of the Forward Calculations 

The basic procedure of our forward technique is to calculate the Stokes vector 
produced along a given LOS in a magnetic model, and build an image from 
a grid of sight lines. To do this, the magnetic field, temperature, density, and 
velocity at every location along each LOS are used. Given this information, 
we calculate the level populations and the emitted polarization profiles for the 
Fe XIII 1074.7 nm transition at each location using a publicly available For- 
tran code (FORCOMP) discussed by Judge and Casini (2001). The forward 
model has an IDL user interface and is publicly available for download and 
use (people.hao.ucar.edu/sgibson/FORWARD/). 

FORCOMP first calculates the statistical-equilibrium equations based on the 
location and the local plasma parameters from the model: height above the solar 
surface [/i], density [p], temperature [T], magnetic field [B], and velocity [v]. 
Using standard atomic data, the statistical equilibrium equations determine the 
relevant level populations of the atomic system for the transition in question. The 
code treats inelastic and superelastic collisional processes, but neglects elastic 
collisions. This omission affects the magnitude, but not the direction, of L and 
leads to a small uncertainty in V. The LOS field strength and the POS field direc- 
tion are not strongly affected by the elastic collisions. Once the level populations 
are determined, FORCOMP solves the radiative-transfer equations to calculate 
the polarization of the reemitted radiation in the direction of the observer (Judge 
and Casini, 2001). The signals are then integrated over wavelength into a single 
number for each pixel (or LOS), and are assembled into an image. 

The benefit of the forward technique is that we can easily calculate the sim- 
ulated polarization signals from a theoretical model of a magnetic system and 
then compare these images with observations. It allows us to test the theories 
against an observable that is directly sensitive to the magnetic field in the corona. 
Additionally it allows for comparison between the models themselves. 

The forward model outputs Stokes /, Q, [/, V and combinations thereof. 
We use mainly intensity [/], relative linear polarization [L/I]^ azimuth [^], and 
relative circular polarization [V/I]. 



4. MHD Models 

For our study of magnetic flux rope and sheared-arcade signatures in the corona, 
we used three models, each having a distinct magnetic morphology: The first 
model is a 3D analytic spheromak flux rope in exact equilibrium (Gibson and 
Low, 1998; Gibson and Low, 2000). The second is an azimuthally symmetric 
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(2.5D) cylindrical flux rope taken from an MHD simulation created to study 
current-sheet formation during CME initiation (Fan and Gibson, 2006). The 
last is a 2.5D sheared arcade taken from MHD simulations of CME initiation by 
the multipolar breakout mechanism (Antiochos, DeVore, and Klimchuk, 1999; 
Karpen, Antiochos, and DeVore, 2012, and references therein). All of the models 
are in, or near, equilibrium and have been argued as models for prominence 
magnetic structure. The two flux-rope models contain a region of concave- up 
magnetic dips that can support prominence plasma against gravity, and they also 
capture many observed properties of coronal cavities (Hudson et al.^ 1999; Gibson 
and Low, 2000; Mackay et al.^ 2010; Reeves et al.^ 2012). Sheared arcades in 2.5D 
typically contain only concave-down regions that can support time-dependent 
prominence condensations if the field is sufficiently flat (Karpen et al.^ 2001). In 
three dimensions, sheared arcades can develop regions of concave- up field lines 
like those in flux ropes, and they then can support the plasma statically against 
gravity (Antiochos, Dahlburg, and Khmchuk, 1994; Luna, Karpen, and DeVore, 
2012, and references therein). 

In the work presented here, we use a single snapshot from the two time- 
dependent MHD models. Because we are studying the steady-state pre-CME 
magnetic structure, all velocities have been set to zero. The times used are 
those where the field is near equilibrium, and thus have close to zero velocity 
everywhere. 

Two of the models that we study are 2.5D. Azimuthal symmetry creates 
structures that are elongated along the LOS. When structures are highly 3D, 
the magnetic information can become smeared along the LOS, making magnetic 
signatures more difficult to identify. We use 2.5D models because the observa- 
tional signatures of the magnetic field we are studying are clear coronal cavities. 
When a coronal filament channel and associated neutral line are along the LOS 
- nearly parallel with the solar equator - a cavity commonly becomes visible 
implying that cavities are elongated along the LOS (Gibson et al.^ 2010). The 
2.5D assumption is thus justified by cavity observations. 

For the spheromak model, we used the density and temperature provided 
by the analytic model. The parameters were chosen such that the density and 
temperature vary only slightly in the calculation domain (Dove et al.^ 2011). 
For the two MHD models, a range of plasma parameters was explored. The 
goal of this work is to study the impact of the magnetic morphology on the 
polarization signatures. We looked at our models both with the original plasma 
distributions from the MHD simulations and with simple spherically symmetric 
plasma profiles. These new plasma distributions are not strictly in equilibrium 
with the magnetic fields. However, they serve the useful purpose of providing a 
means of disentangling those features in the polarization data that are due to 
the magnetic morphology from those that are heavily influenced by the plasma 
parameters. In addition, any effects of changing the plasma distributions on the 
magnetic field structures would be very small, since the prominence and cavity 
are in the low-/3 regime and are nearly in force-free equilibrium in both MHD 
examples. 

Coronal Stokes vectors are calculated for each theoretical system using the 
forward code described in Section 3. We compare these polarization signatures 
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Figure 1. Field-line traces in the analytic spheromak flux-rope model seen from two different 
view points. Prominence material would theoretically sit below the thick black line, which 
traces the flux-rope axis, in the dips of the magnetic fleld. In the forward calculations, we used 
the orientation in b), with the magnetic fleld along the axis pointed toward the observer. An 
animated version of this flgure is available in the electronic supplementary material. 



with each other to determine their similarities and differences and to identify 
their distinguishing features. 

4.1. Model Descriptions 

The first magnetic system we explored is that of an analytic spheromak flux rope. 
A more detailed study of forward model results from this particular flux rope 
can be found in Dove et al. (2011). The spheromak model (Figure 1) is an exact 
solution to the MHD equations in full magnetostatic equilibrium (Gibson and 
Low, 1998). The magnetic field of the flux rope is a closed, twisted- flux system, 
attached to the photosphere, which has been shown to reproduce observational 
features of a three-part CME including the cavity and the bright prominence 
(Gibson and Low, 2000). The external field has a split bipolar configuration 
with a hydrostatic density background. We used an orientation such that the 
flux-rope axis, and hence the prominence material, is oriented along the LOS, 
and the axial magnetic field is directed toward the observer (Figure 1(b)). As 
stated by Dove et al. (2011), we chose a parameter set such that the density is 
close to spherically symmetric. The background-density profile was taken from 
Schmit and Gibson (2011). The density decreases from around 5 x 10^ cm~^ at 
photosphere to about 3 x 10^ cm~^ near the top of the spheromak at 1.3 R©. 
The temperature is between 7 x 10^ and 1 x 10^ K. The magnetic-field strength 
is strongest at the axis where it is around 1 G, and the external field strength 
near the flux rope is of order 0.1 G. Thus, the plasma /3 is high outside the 
spheromak, above 100, and between 1 and 10 in most of the flux rope. 

The second model is a 2.5D axisymmetric cylindrical flux rope (Figure 2) 
and is described in Fan and Gibson (2006). This numerical model comprises a 
2 X 10^ K isothermal atmosphere occupied by a potential arcade field under which 
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Figure 2. Field lines traced from the 2.5D axisymmetric cylindrical flux rope. Colors on the 
fleld lines represent the strength of the y, or axial, component of the fleld. An animated version 
of this flgure is available in the electronic supplementary material. 



a twisted toroidal flux tube is kinematically emerged. With continual emergence, 
the system is quasi-static until time t = 114 Rq/'^ao where '^ao is a characteristic 
Alfven speed, after which the flux rope erupts. If the emergence is stopped before 
that, the system remains stable. We analyzed a time step at t = 114 where the 
emergence was stopped at t = 112. The model extends from 1 R© to 14.4 R© 
radially and from 7r/3 to 27r/3 in latitude. For the forward analysis, the system 
is oriented such that flux rope is in the equatorial plane and the axial field points 
away from the observer. The magnetic-field strength at the axis is around 10 G. 
The electron density in the flux rope is on the order of 10^~^ cm~^. The plasma 
j3 is below 0.1 except near the footprint of the arcade immediately surrounding 
the flux rope. 

The third model (Figure 3) is a 2.5D axisymmetric numerical datacube of 
a breakout quadrupolar system (Karpen, Antiochos, and DeVore, 2012). The 
computational domain extends from 1 to 125 R© radially and over tt radians 
in latitude. The system is energized by shearing the innermost polarities near 
the neutral line. The shearing results in a field that is pointed away from the 
observer in our orientation. The field that connects these two polarities shears 
and expands both radially and laterally. Because a simple adiabatic energy equa- 
tion and a closed lower boundary condition are used, the plasma entrained in 
the sheared field rarefies and adiabatically cools relative to the background. This 
configuration does not form a flux rope with a central axis until flare reconnection 
sets in following the eruption. The time shown in Figure 3 is t = 60 000 seconds; 
stopping the shearing motions at this time results in a stable equilibrium state, 
while continuing the motions leads inexorably to an eruption. In and around the 
sheared arcade, the magnetic-field strength is around 1 G. In this same area. 
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Figure 3. Field lines traced from the 2.5D sheared arcade model. Colors on the field lines 
represent the strength of the y—, or axial, component of the field. An animated version of this 
figure is available in the electronic supplementary material. 



the plasma densities are 10 cm~'^. The temperature reaches a local minimum 
around 3 x 10^ K within the sheared field, while it is about 10^ in the surrounding 
unsheared field. The plasma /3 is below 0.01 in the sheared region, and is of order 
unity in the unsheared region (except of course, in the immediate vicinity of the 
null point, where the /3 becomes very large). 



5. Results 

Interpretation of a polarization image is not necessarily straightforward (see 
Judge, 2007; Rachmeler, Casini, and Gibson, 2012 for detailed descriptions of 
the signal interpretation). For L/I images {i.e. Figure 4(b)), the magnitude of 
the signal is always positive, and the images are usually plotted on a logarithmic 
scale. Bright areas indicate magnetic field that is primarily in the POS. Dark 
areas generally indicate magnetic field that is along the LOS. Sharp, elongated 
dark structures are usually indicative of Van Vleck inversions, marking where 
the magnetic field is at an angle of ^ 54° from radial. The L/I images may have 
magnetic-field direction overlaid as arrows or lines. In the images presented in 
this article, the red arrows indicate the true POS direction of the magnetic field 
in a thin POS slice that bisects the Sun. The blue lines indicate the azimuth 
direction of the linear polarization, which is subject to the 90° Van Vleck am- 
biguity (Section 2.1). Notice that they are generally parallel when the magnetic 
field is close to radial, and perpendicular when it is not. V/I images are plotted 
on a linear scale with blue as negative and red as positive; white is zero. In our 
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coordinate system, positive Stokes V (blue) indicates magnetic field toward the 
observer. 

5.1. Spheromak Flux Rope 

We use the parameter set from Dove et al. (2011) to demonstrate the main 
features for the spheromak flux rope. In the LOS integrated images, the following 
features are identifiable and are robust signatures of this magnetic morphology. 
Figure 4 shows the forward model results, as presented also in Dove et al. (2011). 
We summarize the conclusions from that analysis as follows: 

i) Dark L/I central core. A dark core is clearly visible at the location of the 
flux-rope axis in the L/I image (Figure 4(b); Z = 0, Y = 1.1). This is due 
to the LOS field associated with the axis. The axis of this type of flux rope 
is particularly clear because it is straight along the LOS, and not curved like 
the axis of the cylindrical flux rope. The dark central core is visible if the axis 
is oriented within about 30° of the LOS. 

ii) Dark L/I outer ring. A ring of darker L/I is visible at the edge of the 
spheromak bubble (Figure 4(b)), and it is also associated with LOS field. 
This ring is much fainter than the axis field as the B on the outer edge of the 
bubble is only aligned with the LOS in a relatively narrow volume of space. 

iii) Bright L/I ring. Between i) and ii) is a bright ring in L/I (Figure 4(b)), 
which is due to the POS field in the flux rope. 

iv) Radial azimuth.The linear polarization direction shows no clear Van Vleck 
inversion locations (blue lines in Figure 4(c)). Although there are Van Vleck 
inversions within the spheromak, the rotation of the field along the LOS 
smears these out such that they are not visible in the integration. 

v) Bi-directional circular polarization (V/I). The circular polarization comprises 
a clear circular positive signal around the axis surrounded by a weaker rung 
of negative signal (Figure 4(d)). The presence of both positive and negative 
Stokes V is not found in either of the other models studied here. 

5.2. Cylindrical Flux Rope 

Figure 5 shows the forward-model results for the cylindrical flux-rope config- 
uration. To test the robustness of the magnetic signatures in the polarization 
signals, we ran the forward calculations on three cases. Case Fl uses the original 
density and isothermal temperature of T = 2 x 10^ K from the MHD model 
(Figure 5). Case F2 uses an isothermal hydrostatic density fall-off with a scale 
height of {2RQk})T) / (GMqUIp) « O.ISR© (R© and M© are the solar radius and 
mass, ki) is Boltzman's constant, G is the universal gravity constant, rrip is the 
mass of a proton, and T = 1.5 x 10^ K is the temperature) and a density at the 
coronal base of 5.8 x 10^ electrons cm""^. Case F3 uses a hydrostatic power-law 
density and temperature function derived from fits to coronal-streamer densities 
(Gibson et al.^ 1999). The magnetic field in this model is in nearly force- free 
equilibrium. The imposed plasma is also low-/3. If re-relaxed to true equilibrium, 
the scale-height of the plasma along field lines would be altered, but the field 
topology would remain nearly unchanged. 
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Spheromak: Stokes I (log scale) 



Spheromak: L/l (log scale) 
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Figure 4. Forward-model results of the spheromak flux rope, a) Stokes intensity, b) relative 
linear polarization, c) relative linear polarization with magnetic-fleld direction plotted as red 
arrows and integrated polarization azimuth direction plotted as blue lines, d) relative circular 
polarization. 



Stokes / changed noticeably when a sphericahy symmetric density was used. 
Almost none of the structure seen in Fl (Figure 5(a)) is present in F2 or F3 (not 
shown). This is not surprising as the intensity of emitted radiation is strongly 
dependent on the local plasma density, so a spherically symmetric density results 
in a virtually spherically symmetric Stokes /. The relative linear polarization 
signals had little variation between the three cases (Figure 6); these differences 
are discussed at the end of this Section. We will first analyze the signatures of 
Fl. 

i) Bark V-shaped Van Vleck inversions in the arcade. These are the Van Vleck 
inversions in the external arcade field (Figure 5(b); Y = 1.4 to 2). The field 
outside of the V is less than 54° from radial, and the field inside the V is 
greater than that. This is clearly visible (Figure 5(c)) in the shift of the 
linear-polarization direction (blue lines) from parallel to perpendicular to the 
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Figure 5. Forward-model results from the flux rope model case Fl with the density of the 
original MHD data cube and a temperature of 2 x 10^ K. a) Stokes intensity, b) relative 
linear polarization, c) relative linear polarization with magnetic-fleld direction plotted as red 
arrows and integrated polarization azimuth direction plotted as blue lines, d) relative circular 
polarization. 



Case Fl : L/l (log scale) 



Case F2: L/l (log scale) 



Case F3: L/l (log scale) 




Figure 6. Comparison of the relative linear polarization for the three cases of the cylindrical 
flux rope. All three images use the same scale. 



POS magnetic-field direction (red arrows). These two Van Vleck inversions 
extend downward into the flux rope. 

ii) Darker central region in L/L In general, the region near the flux-rope axis 
(Z = 0, y = 1.35) is slightly darker (Figure 5(b)). This is because there 
is more LOS fleld in this region, so there is less linear-polarization signal. 
However, because of the LOS integration and the curvature of the central 
flux-rope axis, the location of the axis is unclear. See also iii) below. 
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iii) Dark beads in L/I near the axis. These are visible in Fl and F3, but not in F2 
(Figure 6). These locations are dark because they are signatures of the LOS 
field in the flux rope axis. This flux rope is axisymmetric, so it curves around 
the limb of the Sun. The curvature results in a perspective effect such that the 
location of LOS field is only co-spatial with the location of the flux rope axis 
in the central POS slice. The true location of the flux-rope axis is not readily 
apparent in the integrated data because the volume of space where they are 
co-located is a small. The location of the LOS field follows a D-shaped arc 
whose legs coincide with iv). 

iv) Two dark horizontal lines in L/I intersecting the photosphere. These are a 
combination of a Van Vleck inversion in the lower part of the flux rope and 
LOS field from the legs of the D-shape described in iii) and are thus the 
darkest features in the image. 

v) Slightly spoked azimuth direction. The azimuth directions are mostly radial, 
but there is a slight spoke-like signature around the flux-rope axis (Fig- 
ure 5(c)). Since the Van Vleck inversion locations are obvious (in this model), 
the 90° ambiguity can be removed, and the flux rope nature of the field 
becomes evident. Even without removing the ambiguity, the slight spoke may 
be a feature that can help to identify this type of magnetic morphology. Note 
that the linear-polarization azimuth direction (blue lines) everywhere is close 
to radial. 

vi) Bulb of circular polarization. The circular polarization is all the same sign 
(Figure 5(d)). The strongest signal comes from above the hmb, and surrounds 
the location of the flux rope axis. 

Although almost all of the listed features are present in each case, some are more 
pronounced in certain cases. For instance, the dark beads in iii) are distinctly 
visible in F3, and not at all in F2. The density differences in each case change 
the weighting of the signal along the LOS. Thus, certain features are more or 
less clear depending on whether the signal is concentrated at the central POS, 
or spread out along the LOS. Case F2 has the most gradual density drop with 
height and thus the dark beads from iii) are overcome by brighter signal in the 
foreground and background. 

5.3. Sheared Arcade 

Much like the cylindrical flux-rope model, we ran the forward calculations on 
the sheared arcade model with several different density and temperature pro- 
files. The cases presented here are as follows: SO - the density and temperature 
provided by the MHD model; S2 - the density provided by the MHD model and 
isothermal temperature of 1.5 x 10^ K; S3 - hydrostatic streamer density and 
temperature fit from Gibson et al. (1999). We do not relax the configuration 
to equilibrium with the imposed plasma parameters. In all cases, the plasma 
is low-/3 except in the region of the null-line. We will not discuss results from 
SO. The temperature of the original data falls below the minimum temperature 
threshold for the forward calculations (around 5 x 10^ K). Most of the sheared- 
field plasma is at or below this temperature, due to the assumption of adiabatic 
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energy transport, so any calculations on the original data only produce signal 
from the unsheared field, which is not useful for this work. 

Figure 7(g), (h) shows the comparison of the integrated L/I for the S2 and 
S3 cases in the sheared region. We have looked at the polarization signals from 
thin POS slices along the LOS, and find that for any given slice, the L/I is 
qualitatively similar between the two cases. The differences seen in Figure 7 
arise from the relative contributions to the integrated signal from the sheared 
field versus the background unsheared field. Stokes /, Q, and U are weighted by 
density. Hence, higher-density regions contribute more to the integrated signal 
than lower density regions. In S2, the original density profile is used, and in S3 
a spherically symmetric density is used. These densities are plotted in Figure 
7(a), (b). Note that in S2, the density in the unsheared field is about an order 
of magnitude more dense than the sheared field due to the large volumetric 
expansion of the sheared field and the closed lower boundary condition imposed 
in the simulation. From a broad range of observations, it is known that streamers 
tend to be about a factor of two more dense than the embedded cavity (Fuller 
and Gibson, 2009; Schmit and Gibson, 2011). Thus, the S3 results put too much 
emphasis on the sheared field, and the S2 results put too little emphasis on this 
region compared to observations. 

The following list describes the features that are present in the data. Where 
certain features are found to a lesser degree, or not found at all, in one of the 
cases, it is explained in more detail after the list. 

i) Quadrupolar Van Vleck signal. A clear Van Vleck signature is associated with 
the quadrupolar field (Figure 7(c), (d)). Even in the absence of any sheared 
field, there would be three pairs of Van Vleck inversions associated with the 
inner three loop systems. These pairs are the top two, the middle two, and 
the lower two elongated Van Vleck nulls. 

ii) V-shaped Van Vleck inversion. At the top of the central loop system, there are 
Van Vleck lines in a V-shape (Figure 7(e), (f); at a height of 1.6 — 2 R©), which 
is analogous to property i) in the cylindrical fiux-rope model. The sheared 
field is confined to the central part of this system {i.e. the areas of negative 
V/I in Figure 7(i),(j)) and the V-shaped Van Vleck lines are associated with 
the unsheared portion of this central magnetic-loop system (dark-blue loops 
above the sheared field in Figure 3). 

iii) Parallel Van Vleck inversions. In the sheared field, there are two dark parallel 
Van Vleck inversion lines (Figure 7(e), (f); at a height of 1 — 1.6 R©) that 
emanate from the photosphere and connect to the V-shape listed in ii). The 
parallel Van Vleck inversion lines are associated with the legs of the sheared 
region. This field is inclined toward/away from the observer at ^ 54° from 
radial. These inversions are more pronounced in S3 than in S2. 

iv) Dark LOS core in L/L The central sheared area has generally lower linear 
polarization magnitude due to the presence of fields that are more LOS than 
the surrounding field. This effect is stronger in S3 than in S2. 

v) Anomalous LOS signal in L/L A dark spot in L/I is visible in the sheared 
field region in case S3 (Figure 7(f); Z = 0, F = 1.35). This is not due to an 
axis of LOS field. This anomalous LOS signal (Rachmeler, Casini, and Gibson, 
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Figure 7. Forward-model results from the sheared arcade system. The left column shows case 
S2 and the right column shows case S3, a) and b) density in the central region, the electron 
number density is given by 10^ cm""^ where X is the value indicated by the color bar; c) 
and d) relative linear polarization of the entire system; e) and f ) relative linear polarization of 
the central region; g) and h) same as above overlaid with red field vectors and blue azimuth 
directions; i) and j) relative circular polarization of the inner region. 
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2012) arises from cancellation in the LOS integrated Stokes Q and U due to 
the symmetry in the system and is dependent on the relative density in the 
sheared and unsheared regions. This feature is not present in S2 (Figure 7(e)). 

vi) Non-radial azimuth direction. In the sheared-field region, the linear polar- 
ization azimuth direction is parallel to the limb in the S3 case (blue lines 
in Figure 7(h)), which creates an illusion of a cylindrical flux rope. This is 
an area where the magnetic-field angle exceeds the Van Vleck angle, so the 
true field is perpendicular to the observed azimuth direction. Note that at 
this location, the field is actually predominantly in the LOS, and the POS 
component is small. This feature is not present in S2 (Figure 7(g)). 

vii) Strongest circular polarization near limb (V/I). The relative circular polariza- 
tion is shown in Figure 7(i),(j). The shear is not concentrated above the limb, 
as in the flux-rope case, but at the limb. The features within the negative 
V/I in S2 are due to density variations (Figure 7(a)). The V/I in the S3 case 
shows a smooth profile. 

Most of these features are present in both cases, but they are not always obvious 
in the integrated S2 data. The locations of the Van Vleck inversion and the 
LOS field i)-iii) are the same in both cases. Figure 8 shows the LOS integrated 
values of Stokes Q and U [L = \JQ^ + 'U^\ for cases S2 and S3. The locations of 
lowest Q and U are the same for both cases, but in S2, there are no true nulls in 
Q inside the sheared-field region. Thus, the integrated L/I for S2 (Figure 7(e)) 
does not show true Van Vleck nulls. The inversions are there, but the signal from 
the unsheared portion of the LOS obscures them. 

There is clearly no anomalous LOS signal iv) in S2 (Figure 7(e)). This is 
because the background dominates the integrated signal at all heights in this 
case. Investigation of the polarization generated in thin POS slices shows that the 
unsheared field in all cases produces a negative Stokes Q signal, and the sheared 
field produces signal that is predominantly positive. The relative density in the 
sheared and background field dictates where the sheared field dominates in the 
integrated Stokes profiles. The integrated Stokes Q signal for S3 (Figure 8(b)) 
is positive in the central region, showing that the sheared field dominates the 
LOS signal there. Where a zero-line in Stokes Q crosses a zero-line in Stokes U 
(Figure 8(b),(d)) there is an anomalous LOS signal (Figure 7(f)) (Rachmeler, 
Casini, and Gibson, 2012). Since the sheared field in S2 has less density, the 
integrated Stokes Q is always dominated by the background field (Figure 8(a)) 
and is never negative. Thus there is no anomalous LOS signal in S2. 

The relative weighting of the sheared versus the unsheared field also causes the 
difference in the azimuth direction for cases S2 and S3. The volume of space that 
contains sheared field lies inside a Van Vleck inversion - such that the azimuth 
is perpendicular to the POS magnetic-field direction - and inspections of a thin 
POS slice reveals azimuth directions that are consistent with Figure 7(h). The 
integrated L/I in case S2 shows radial azimuths because the background field 
dominates, so the non-radial azimuth signal is overwhelmed. 
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Figure 8. Integrated Stokes Q and U signals from cases S2 and S3 of the sheared arcade 
model. Lines of zero Q or U are shown in magenta. 



6. Discussion 

We have presented three coronal models and their synthetic polarization signa- 
tures. We find that each of these models is distinct and distinguishable, even 
when using linear polarization alone. The spheromak fiux rope is the most 
recognizably different, while the cylindrical fiux rope and the sheared arcade 
models have some similarity. 

This work highlights the importance of using a forward approach on coronal 
emission-line polarization. It can teach us what to look for in observations, such 
as Van Vleck inversions and LOS field. It also calls attention to the fact that 
we cannot trust our intuition to pick out magnetic morphologies. The sheared- 
arcade model is a good example of this. On initial inspection, the polarization 
signatures of the S3 sheared arcade resemble a fiux rope; the azimuth direction 
is parallel to the limb between the inner Van Vleck inversions, and a false axial 
signature may be present. Both of these signatures can be logically explained 
when analyzing the forward results, but this example underlines the need for 
forward or inverse analysis before magnetic- structure identification can be made. 
Another important strength of the forward approach is that it fully takes into 
account the LOS integration of the polarization signal. The presence of the 
optically thin plasma is a significant challenge for the inverse technique and 
so is often seen as a limitation for coronal polarization data as a whole, but 
the forward approach incorporates the lack of opacity. By looking at a given 
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magnetic configuration with multiple plasma profiles, we can also learn about 
how the signatures change due to the plasma parameters alone. We have shown 
that in the case of cavities, where the cavity itself is about half as dense as the 
surrounding streamer, the polarization signature from the streamer can obscure 
some of the features that are present in the cavity. For future observations, it 
is clear that knowledge of the LOS density structures is important for analysis 
and interpretation. Density diagnostics which determine a 3D density distribu- 
tion could be used in conjunction with polarization observations and forward 
calculations. 

Our work is not the first to use this forward approach to understand hypo- 
thetical or actual observations. Judge, Low, and Casini (2006) applied the same 
forward code to study prominence-supporting magnetic fields and current sheets. 
Liu and Lin (2008) compared observations of an active region on the limb to 
potential-field extrapolations to study how the LOS affects the fit of the forward 
calculation to the observations. Dove et al. (2011) compared the spheromak 
model presented here to an early CoMP observation of a large cavity. The next 
important step is to take the knowledge gained with these forward studies and 
apply it directly to observations, looking for the specific morphologies; this work 
is already underway. Bak-Steshcka et al. (2013) have found that cavities observed 
in CoMP in 2011 and 2012 usually have a characteristic "rabbit-head" signature 
in L/I. This signature consists of two Van Vleck inversions, with or without a 
dark central region indicating LOS field. They have shown that this observation 
is consistent with a 3D flux-rope topology where the height of the dark central 
"head" is approximately co-spatial with the center of the cavity. 

Observations carry their own challenges because of noise and the existence of 
small-scale density structures in the corona. CoMP is an occulted coronagraph 
and the occulter is at approximately 1.05 R©, which means that especially 
for small cavities, the distinguishing characteristics such as azimuth direction 
would likely be obscured by the occulter. We have used extremely simplified 
density structures to isolate the magnetic features, but in reality, the Stokes 
/ observations are highly structured. By analyzing relative linear and circular 
polarization, we remove some of the density component, but the signal is still 
density dependent and we have shown that the relative importance of the struc- 
tures along the LOS is highly dependent on their density. Our current approach 
is more applicable in coronal cavities which are, in general, fairly smooth in 
intensity compared to active regions with clear bright loops. In future forward- 
modeling research, more realistic density models are needed. The observational 
noise, the occulter, and the highly structured coronal density make it difficult 
to uniquely characterize observed cavities, as they may be consistent with more 
than one model. 

We are just beginning to scratch the surface of what the polarization data 
can teach us about the solar corona. Here we have studied idealized equilibrium 
structures. Not only is there a large range of magnetic morphologies left to 
study, there is also the important aspect of time-dependence that is still open 
for exploration. The forward approach is only one of the methods available for 
analyzing these data, and there is still much to do with comparisons to obser- 
vations, true forward fitting for given observations, and looking at the Sun as a 
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whole as opposed to specific magnetic structures. We look forward to witnessing 
the advances that come out of these data in conjunction with both forward and 
inverse techniques. 
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